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Composition of Low Redshift Halo Gas 

Renyue Cen^ 

ABSTRACT 

Halo gas in low-z (z < 0.5) > O.IL^, galaxies in high- resolution, large-scale 
cn cosmological hydrodynamic simulations is examined with respect to three com- 

Q ponents: (cold, warm, hot) with temperatures equal to (< 10^, 10^~^, > 10^)K, 

CN respectively. The warm component is compared, utilizing O VI AA1032, 1038 

Q_^ absorption lines, to observations and agreement is found with respect to the 

-^ galaxy-0 VI line correlation, the ratio of O VI line incidence rate in blue to 

T-H red galaxies and the amount of O VI mass in star-forming galaxies. A detailed 

account of the sources of warm halo gas (stellar feedback heating, gravitational 
shock heating and accretion from the intergalactic medium), inflowing and out- 
flowing warm halo gas metallicity disparities and their dependencies on galaxy 
types and environment is also presented. Having the warm component securely 
p I anchored, our simulations make the following additional predictions. First, cold 

Q gas is the primary component in inner regions, with its mass comprising 50% of 

^— • all gas within galacto-centric radius r = (30, 150)kpc in (red, blue) galaxies. Sec- 

c3 ond, at r > (30, 200)kpc in (red, blue) galaxies the hot component becomes the 

majority. Third, the warm component is a perpetual minority, with its contribu- 
^ tion peaking at ~ 30% at r = 100 — 300kpc in blue galaxies and never exceeding 

^ 5% in red galaxies. The significant amount of cold gas in low-z early-type galax- 

■^ ies found in simulations, in agreement with recent observations (Thom et al.) is 

^ intriguing, so is the dominance of hot gas at large radii in blue galaxies. 

O 

cn Subject headings: hydrodynamics — shock waves — supernovae: general — 

galaxies: halos — intergalactic medium — cosmology: theory 

X 

1. Introduction 

Galaxy formation and evolution is the central astrophysical problem in cosmology. 
The basic parameters of the cosmological framework - the standard cosmological constant- 



dominated cold dark matter (DM) model (LCDM) (e.g., Krauss & Turner 1995 Bahcall 
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et al. 1999) - are largely fixed to an accuracy of ~ 10% or better. The LCDM model is able 



to explain a variety of observations on scales greater than ~ IMpc, including high redshift 



supernovae (e.g., Perlmutter et al. 1998 Riess et al. 1998 Astier et al. 2006), the cosmic 



microwave background (e.g., Komatsu et al. 2011 Planck Collaboration et al. 2013), large 



scale distribution of galaxies (e.g., Tegmark et al. 2004 Percival et al. 2007), X-ray cluster 
abundance (e.g., [Allen et al.||2008D and Lya forest (e.g., [Croft et al.||2002[[Seljak et al.||2005D . 



An important component of the astrophysical problem - gravitational formation and 
evolution of halos that host galaxies - is well understood, through N-body simulations (e.g., 
Jenkins et al] [20011 [Bullock et al]|2001t [Wechsler et~al]|2002[ [Diemand et al.[[2007D and 
analytic models (e.g., [Bond et al.|[1991| [Lacey & Cole|[1993| [Sheth & Tormen[[1999| [Mo & 



White 2002 Cooray & Sheth 2002). The gastrophysics of galaxy formation and feedback. 



on the other hand, is far from being adequately understood. Alternative approaches that 
parameterize and then infer physical processes based on finding best matches to observations, 
such as the semi-analytic methods (e.g., Somerville fc Primack[[l999 Benson et al.[[2003 ) and 
the halo-occupation distribution (HOD) method (e.g., Berlind fc Weinberg|2002 Zheng et al. 



2007), have been successful but have limited predictive power. More importantly, in semi- 



analytic methods the treatment of galaxy formation is halo based and largely decoupled 
from that of the intergalactic medium, which in fact has dramatically evolved with time. 
At z = 2 — 6 most of the baryons are found to be in the Lya forest, a relatively cold 
phase of temperature of ~ lO'^K, as indicated by both observations (e.g.. Ranch et al. 1997) 
and simulations (e.g., Cen et al. 1994). By z = most of the baryons in the intergalactic 



medium have been heated up, primarily by gravitational shocks, to temperatures that are 
broadly peaked at about lO^K, the so-called Warm- Hot Intergalactic Medium (WHIM) (e.g., 
Cen & Ostriker 1999). The "ab initio", more predictive approach of direct cosmological 



hydrodynamic simulations, after having made steady progress (e.g., Evrard et al.[[1994 Katz 
eTaLlpigGl [Teyssier[[2002| [Keres et al.|[2005| [Hopkins et aT][2006l [Oppenheimer fc Dave[[2006 



Governato et al. 2007 Naab et al. 2007 Gnedin et al. 2008 Joung et al. 2009 Cen 2011) 



begin to be able to make statistically significant and physically realistic characterizations of 
the simultaneous evolution of galaxies and the intergalactic medium. 

It is the aim of this writing to quantify the composition of the halo gas in low redshift 
galaxies, using state-of-the-art high resolution (460/i~^pc), large-scale (thousands of galax- 
ies) cosmological hydrodynamic simulations with advanced treatments of star formation, 
feedback and microphysics. Our focus here is on gas that is in the immediate vicinities of 
galaxies, on galactocentric distances of 10 — SOOkpc, where the exchanges of gas, metals, 
energy and momentum between galaxies and the intergalactic medium (IGM) primarily take 
place. We shall broadly term it "circumgalactic medium (CGM)" or "halo gas". Under- 
standing halo gas is necessary before a satisfactory theory of galaxy formation and evolution 
may be constructed. The present theoretical study is also strongly motivated observation- 
ally, in light of recent rapid accumulation of data by HST observations enabling detailed 



comparisons between galaxies and the warm component (T ~ 10^ — lO^K) of their CGM at 



low redshift (e.g., Chen & Mulchaey 2009 Prochaska et al. 2011b; Tumlinson et al. 2011a 



Tripp et al. 2011). 



We shall dissect halo gas at low redshift {z < 0.5) into three components, (cold, warm, 
hot) gas with temperature (< 10^, 10^ — 10®, > 10®)K, respectively. A large portion of 
our presentation is spent on quantifying O VI AA1032, 1038 absorption hues and the overall 
properties of warm halo gas and comparing them to observations in as much detail as possible. 
Feedback processes, while being treated with increased physical sophistication, are still not 
based on first principles due primarily to resolution limitations in large-scale cosmological 
simulations. Thus, it is imperative that our simulations are well validated and anchored by 
requiring that some key and pertinent aspects of our simulations match relevant observations. 
The O VI line, when coUisionally ionized, has its abundance peaked at a temperature of 
T = 10^'^~^''^K and thus is an excellent proxy for the the warm gas. After validating our 
simulations with respect to the observed properties of O VI absorption lines, we present 
the overall composition of low redshift halo gas. We find that, for (red, blue) galaxies more 
luminous than O.IL* the cold gas of T < lO^K, on average, dominates the halo gas budget 
within a radius of (30, 150)kpc. Beyond a radius of (30, 200)kpc for (rcd,blue) galaxies the 
hot gas of T > lO^K dominates. The warm component remains a smallest minority at all 
radii, peaking at ~ 30% at ~ 100 — 300kpc for blue galaxies but never exceeding 5% for red 
galaxies. 

The following physical picture emerges for the physical nature of the warm gas com- 
ponent. The warm halo gas has a cooling time much shorter than the Hubble time and 
hence is "transient", with their presence requiring sources. To within a factor of two we 
find that, for low-z > O.IL^, red galaxies contributions to warm halo gas from star formation 
feedback (F,.), accretion of intergalactic medium (A.,.) and gravitational shock heating (G,.) 
are (i^,-, A^., G,) = (30%, 30%, 40%). For blue > O.IL* galaxies contributions to warm halo 
gas from the three sources are {Ft,, A^, Gb) = (48%, 48%, 4%). The mean metallicity of warm 
halo gas in (red, blue) galaxies is (~ 0.25 Zq, ~ 0.11 Zq). Environmental dependence of 
O Vl-bearing halo gas is as follows. In low density environments the metallicity of inflowing 
warm gas is substantially lower than that of outflowing warm gas; the opposite is true in 
high density environments. 

The outline of this paper is as follows. In §2.1 we detail simulation parameters and hy- 
drodynamics code, followed by a description of our method of making synthetic O VI spectra 
in §2.2, which is followed by a description of how we average the two separate simulations 
C (cluster) and V (void) run in §2.3. Results are presented in §3. A detailed comparison 
of galaxy-0 VI absorber correlation is computed and shown to match observations in §3.1, 
followed in §3.2 by an analysis of the ratio of O VI absorber incidence rates around blue and 
red galaxies that is found to be consistent with observations. A detailed examination of the 



physical origin and properties of the warm gas in low-z halo is given in §3.3. The overall 
composition of low-z halo gas is given in §3.4 and conclusions are summarized in §4. 



Simulations 



2.1. Hydrocode and Simulation Parameters 



We perform cosmological simulations with the AMR Eulerian hydro code, Enzo (Bryan 



1999 Bryan & Norman 1999 O'Shea et al. 2005). The version we use is a "branch" ver- 



sion (Joung et al. 2009), which includes a multi-tiered refinement method that allows for 



spatially varying maximum refinement levels, when desired. This Enzo version also includes 
metallicity-dependent radiative cooling extended down to 10 K, molecular formation on dust 
grains, photoelectric heating and other features that are different from or not in the public 
version of Enzo code. We use the following cosmological parameters that are consistent with 
the WMAPT-normalized ( [Komatsu et~aL][20lI| LCDM model: ^m = 0.28, n^ = 0.046, 
^A = 0.72, as = 0.82, Hq = 100/ikms"^Mpc"^ = TOkms'^Mpc"^ and n = 0.96. These 



parameters are also consistent with the latest Planck results (Planck Collaboration et al. 



2013), if one adopts the Hubble constant that is the average between Planck value and those 



derived based on SNe la and HST key program ( Riess et"aL]|2011 Freedman et al.||2012 ). We 
use the power spectrum transfer functions for cold dark matter particles and baryons using 



fitting formulae from Eisenstein & Hut (1998). We use the Enzo inits program to generate 
initial conditions. 

First we ran a low resolution simulation with a periodic box of 120 h^^Mpc on a side. 
We identified two regions separately, one centered on a cluster of mass of ~ 2 x 10^^ Mq and 
the other centered on a void region at z = 0. We then resimulate each of the two regions 
separately with high resolution, but embedded in the outer 120/i~^Mpc box to properly take 
into account large-scale tidal field and appropriate boundary conditions at the surface of 
the refined region. We name the simulation centered on the cluster "C" run and the one 
centered on the void "V" run. The refined region for "C" run has a size of 21 x 24 x 
20h~^Mpc^ and that for "V" run is 31 x 31 x 35/i~^Mpc^. At their respective volumes, 
they represent 1.8cr and — l.Oo" fluctuations. The root grid has a size of 128^ with 128^ 
dark matter particles. The initial static grids in the two refined boxes correspond to a 
1024^ grid on the outer box. The initial number of dark matter particles in the two refined 
boxes correspond to 1024^ particles on the outer box. This translates to initial condition in 
the refined region having a mean interparticle-separation of 117/i^^kpc comoving and dark 
matter particle mass of 1.07 x lO®/i^^M0. The refined region is surrounded by two layers 
(each of ~ lh~^Mpc) of buffer zones with particle masses successively larger by a factor of 8 
for each layer, which then connects with the outer root grid that has a dark matter particle 



mass 8^ times that in the refined region. The initial density fiuctuations are included up 
to the Nyquist frequency in the refined region. The surrounding volume outside the refined 
region is also followed hydrodynamically, which is important in order to properly capture 
matter and energy exchanges at the boundaries of the refined region. Because we still can 
not run a very large volume simulation with adequate resolution and physics, we choose 
these two runs of moderate volumes to represent two opposite environments that possibly 
bracket the universal average. 

We choose a varying mesh refinement criterion scheme such that the resolution is always 
better than 460/h proper parsecs within the refined region, corresponding to a maximum 
mesh refinement level of 9 above z = 3, of 10 at z = 1 — 3 and 11 at z = — 1. The simulations 
include a metagalactic UV background ( Haardt fc Madau]|2012 ), and a model for shielding of 
UV radiation by atoms ( Cen et al.|2005 ). The simulations also include metallicity-dependent 
radiative cooling and heating ( Cen et al.|1995 ). We clarify that our group has included metal 



cooling and metal heating (due to photoionization of metals) in all our studies since Cen 
et al. (1995) for the avoidance of doubt (e.g., Wiersma et al.|2009 Tepper-Garcia et al.|2011 ). 
Star particles are created in cells that satisfy a set of criteria for star formation proposed by 



Cen & Ostriker (1992). Each star particle is tagged with its initial mass, creation time, and 
metallicity; star particles typically have masses of ~10^~^ M©. 



Supernova feedback from star formation is modeled following Cen et al. (2005 ). Feedback 



energy and ejected metal-enriched mass are distributed into 27 local gas cells centered at 
the star particle in question, weighted by the specific volume of each cell (i.e., weighting 
is equal to the inverse of density), which is to mimic the physical process of supernova 
blastwave propagation that tends to channel energy, momentum and mass into the least 
dense regions (with the least resistance and cooling). We allow the whole feedback processes 
to be hydrodynamically coupled to surroundings and subject to relevant physical processes, 
such as cooling and heating, as in nature. The extremely inhomogeneous metal enrichment 
process demands that both metals and energy (and momentum) are correctly modeled so 
that they are transported into right directions in a physically sound (albeit still approximate 
at the current resolution) way, at least in a statistical sense. In our simulations metals are 
followed hydrodynamically by solving the metal density continuity equation with sources 
(from star formation feedback) and sinks (due to subsequent star formation). Thus, metal 
mixing and diffusion through advection, turbulence and other hydrodynamic processes are 
properly treated in our simulations. 

The primary advantages of this supernova energy based feedback mechanism are three- 
fold. First, nature does drive winds in this way and energy input is realistic. Second, 
it has only one free parameter esN, namely, the fraction of the rest mass energy of stars 
formed that is deposited as thermal energy on the cell scale at the location of supernovae. 
Third, the processes are treated physically, obeying their respective conservation laws (where 



they apply), allowing transport of metals, mass, energy and momentum to be treated self- 
consistently and taking into account relevant heating/cooling processes at all times. We use 
GsN = 1 X 10~^ in these simulations. The total amount of explosion kinetic energy from 
Type II supernovae with a Chabrier IMF translates to e^Ar = 6.6 x 10~^. Observations 
of local starburst galaxies indicate that nearly all of the star formation produced kinetic 



energy (due to Type II supernovae) is used to power galactic superwinds (e.g., Heckman 
2001). Given the uncertainties on the evolution of IMF with redshift (i.e., possibly more 



top heavy at higher redshift) and the fact that newly discovered prompt Type I supernovae 
contribute a comparable amount of energy compared to Type II supernovae, it seems that 
our adopted value for esN is consistent with observations and physically realistic. The 



validity of this thermal energy-based feedback approach comes empirically. In Cen (2012b) 
the metal distribution in and around galaxies over a wide range of redshift (z = — 5) is 
shown to be in excellent agreement with respect to the properties of observed damped Lya 



systems (Rafelski et al. 2012), whereas in Cen (2012a) we further show that the properties 



of O VI absorption lines at low redshift, including their abundance, Doppler-column density 
distribution, temperature range, metallicity and coincidence between O VII and O VI lines. 



are all in good agreement with observations (Danforth &; ShuU 2008 Tripp et al. 2008 Yao 



et al. 2009). This is non-trivial by any means, because they require that the transport of 



metals and energy from galaxies to star formation sites to megaparsec scale be correctly 
modeled as a function of distance over the entire cosmic timeline, at least in a statistical 
sense. 



2.2. Simulated Galaxy Catalogs 



We identify galaxies in our high resolution simulations using the HOP algorithm (Eisen- 



stein & Hu 1999), operated on the stellar particles, which is tested to be robust and insen- 



sitive to specific choices of concerned parameters within reasonable ranges. Satellites within 
a galaxy are clearly identified separately. The luminosity of each stellar particle at each of 
the Sloan Digital Sky Survey (SDSS) five bands is computed using the GISSEL (Galaxy 



Isochrone Synthesis Spectral Evolution Library) stellar synthesis code (Bruzual & Chariot 



2003), by supplying the formation time, metallicity and stellar mass. Collecting luminosity 



and other quantities of member stellar particles, gas cells and dark matter particles yields 
the following physical parameters for each galaxy: position, velocity, total mass, stellar mass, 
gas mass, mean formation time, mean stellar metallicity, mean gas metallicity, star formation 
rate, luminosities in five SDSS bands (ugriz) and others. At a spatial resolution of proper 
460pc/h with more than 2000 well resolved galaxies at z = 0, this simulated galaxy catalog 
presents an excellent (by far, the best available) tool to study circumgalactic medium around 
galaxies at low reshift. 



In some of the analysis we perform here we divide our simulated galaxy sample into two 
sets according to the galaxy color. We shall call galaxies with g — r < 0.6 blue and those 
with g — r > 0.6 red. It is found that g — r = 0.6 is at the trough of the galaxy bimodal 



color distribution of our simulated galaxies (Cen 2011 Tonnesen & Cen 2012|, which agrees 



well with that of observed low-z galaxies (e.g., Blanton et al. 2003). 



2.3. Generation of Synthetic O VI Absorbers 

The photoionization code CLOUDY (Ferland et al. 1998) is used post-simulation to 
compute the abundance of O VI, adopting the shape of the UV background calculated by 



Haardt & Madau (2012) normalized by the intensity at 1 Ryd determined by ShuU et al. 



(1999) and assuming ionization equilibrium. We generate synthetic absorption spectra given 



the density, temperature, metallicity and velocity fields in simulations. Each absorption line 
is identified by the velocity (or wavelength) interval between one downward-crossing and 
the next upward-crossing points at flux equal to 0.99 (flux equal to unity corresponds to an 
unabsorbed continuum flux) in the spectra. We do not add instrumental and other noises 
to the synthetic spectra. Since the absorption lines in question are sparsely distributed in 
velocity space, their identifications have no significant ambiguity. Column density, equivalent 
width, Doppler width, mean column density weighted velocity and physical space locations, 
mean column density weighted temperature, density and metallicity are computed for each 
line. We sample the C and V run, respectively, with 72, 000 and 168, 000 random lines of 
sight at z = 0, with a total pathlength of Az ~ 2000. A total of ~ 30, 000 > 50 mA O VI 
absorbers are identified in the two volumes. While a detailed Voigt profile fitting of the fiux 
spectrum would have enabled closer comparisons with observations, simulations suggest that 
such an exercise does not necessarily provide a more clarifying physical understanding of the 
absorber properties, because bulk velocities are very important and velocity substructures 



within an absorber do not necessarily correspond to separate physical entities (Cen 2012a) 



2.4. Averaging C and V Runs 

The C and V runs at z = are used to obtain an "average" of the universe. This 
cannot be done precisely without much larger simulation volumes, which is presently not 
feasible. Nevertheless, we make the following attempt to obtain an approximate average. 
The number density of galaxies with luminosity greater than O.IL^, in SDSS r-band in the 
two runs is found to be 3.95 x 10~^/i^Mpc"^ and 1.52 x 10"^/i^Mpc"^, respectively, in the C 
and V box. We fix the weighting for C and V run for the purpose of averaging statistics of 
the C and V runs by requiring that the average density of galaxies with luminosity greater 
than O.IL^ in SDSS r-band in the simulations to be equal to the observed global value of 



2.87 X lO-^h-^Mpc"-^ by SDSS (Blanton et al. 2003). In the results shown below we use this 



method to obtain averages of statistics, where doing so allows for some more quantitative 
comparisons with observed data. 



3. Results 



3.1. Galaxy-O VI Absorber Correlation at z = — 0.5 



W>50mA. >0.1 L., z=0-0.2 (full sim; 2a error) 
.W>50mA, >0.1 L., z=0-0.2 (photoioniz; 2a error) 

♦ Chen & Mulchaey (2009): z=0-0.5 (obs) 
D Prochaska et al (2011 ) : z=0-0.2 (obs) 

• Tumlinson etal (2011): z=0.1-0.36 (obs) 
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Fig. 1. — Cumulative probability distribution functions of > 50 mA O VI absorbers of 
finding > O.IL* galaxies at z = — 0.2 from simulations with 2a errorbars (red solid curves). 
The distribution functions at z = — 0.2 are obtained by averaging z = and z = 0.2 results 



with equal weighting. Also shown as symbols are observations from Chen & Mulchaey (2009) 



(solid diamonds), Prochaska et al. (2011a) (open squares) and Tumlinson et al. (2011b) 



(solid dots). Because the impact parameter of Tumlinson et al. (2011b) samples reaches 
only 150kpc, we have normalized their data points by matching their rp = 150kpc point to 
the Tp = 150kpc point of Prochaska et al. (2011a). The blue dashed curve is produced when 



only photoionized O VI lines with temperature T < 3 x lO^K in our simulations are used. 
The X square per degree of freedom for the red solid curve using all observed data points is 
1.2, whereas it is 7.6 for the blue dashed curve. 



Figure [T] shows the cumulative probability distribution functions of > 50 mA O VI 
absorbers of finding > O.IL^, galaxies at z = — 0.2 from simulations as well as observations. 



We find good agreement between simulations and observations, quantified by the x square per 
degree of freedom of 1.2. In comparison, if using only the low temperature (T < 3 x lO'^K) 
O VI absorbers in the simulations, the cumulative probability is no longer in reasonable 
agreement with observations, with the x square per degree of freedom equal to 7.6; this 
exercise, however, only serves as an illustration of what a photoionization dominated model 
may produce. It will be very interesting to make a similar calculation directly using SPH 
simulations that have predicted the dominance of photoionized O VI absorbers even for 
strong O VI absorbers as shown here (e.g.. 



Oppenheimer et al. 2012) 



This significant difference found with respect to the strong O VI absorber-galaxy cross 
correlations between the photoionization and coUisional ionization dominated models stems 
from the relative difference in the locations of strong O VI absorbers in the two models. In 



the coUisional ionization dominated model (Cen 2012a ShuU et al. 2011) the strong O VI 



absorbers are spatially closer to galaxies in order to have high enough temperature (hence 
high O VI abundance) and high enough density to make strong O VI absorbers, whereas in 



the photoionization dominated model (Tepper-Garcia et al. 2011 Oppenheimer et al. 2012) 



they have to be sufficiently far from galaxies to have low enough densities to be photoionized 
to O VI. Additional requirement in the latter for production of strong O VI absorbers is 
high metallicity (> 0.1 Z©) to yield high enough O VI columns, as found in SPH simulations 



(Tepper-Garcia et al. 2011 Oppenheimer et al. 2012) 



3.2. O VI Absorbers Around Blue and Red Galaxies 



Observations have shown an interesting dichotomy of O VI incidence rate around blue 
and red galaxies. Figure |2] shows the cumulative probability distribution functions of A^ovi > 
10^^ cm~^ O VI absorbers of finding a (red, blue) galaxy of luminosity of > O.IL,,, at z = 0.2 
from simulations. In Figure [3] we show the ratio of the cumulative radial distribution per 
red galaxy to per blue galaxy of > O.IL^, at z = 0.2, compared to observations. It is seen in 
Figure [Sjthat in the r = 50 — 300kpc range the ratio of incidence rate of strong O VI absorbers 
around red galaxies to that around blue galaxies is about 1:5, in quantitative agreement 
with observations. In the case with photoionized O VI absorbers only, the fraction of O VI 
absorbers around red galaxies is much lower and lies significantly below the observational 
estimates, although the present small observational sample prevents from reaching strong 
statistical conclusions based on this ratio alone. 

We see in Figure |2] that statistical uncertainties of the radial probability distribution of 
simulated O VI absorbers are already very small due to a significant number of simulated 
galaxies and a still larger number of simulated absorbers used. What limits the ability to 
make firm statistical statements is the sample size of observational data. Hypothetically, if 
the mean remains the same, a factor of two smaller errorbars would render the photoioniza- 
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Fig. 2. — The cumulative probability distribution functions of O VI absorbers with column 
density greater than 10^^ cm~^ of finding a (red, blue) galaxy of luminosity of > O.IL* (in 
SDSS r-band) aX z = 0.2 with (red dashed, blue solid) curves from simulations with lOcr 
errorbars. The (red dotted curve, blue dot-dashed curve) are the corresponding functions 
for the subset of O VI absorbers that have temperature T < 3 x lO^K in our simulations. 

tion dominated model inconsistent with observations at > 2a confidence level, whereas our 
coUisionally dominated model would be consistent with observations within la. 



3.3. Physical Origin of O VI Absorbers 

We now turn to an analysis to give a physical description for the origin of O VI absorbers 
in the CGM, in the context of the cold dark matter based model. While this section is 
interesting on its own for physically understanding halo gas, the next section on halo gas 
mass decomposition is not predicated on it. 

First, we ask whether the warm gas traced by O VI absorbers requires significant 
energy input to be sustained over the Hubble time. Figure |4] shows the metal mass in 
the warm gas (T = 10^ — lO^K) distributed in the density-metallicity phase space. We 
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Fig. 3. — The ratio of cumulative radial probability distribution function of O VI absorbers 
of equivalent width (W) greater than SOmA per red galaxy to that per blue galaxy of > O.IL* 
(in SDSS r-band) at redshift z = 0.2 from simulations (the black solid curve, 2a errorbars). 
The same ratio for photoionized O VI absorbers (T < 3 x lO^K) only is shown as the green 



dashed curve. Also shown as an open circle is the observation by Tumlinson et al. (2011b). 



see that most of warm metals is concentrated in a small phase space region centered at 
(n, Z) = (10~^cm~^, 0.15 Zq). We note that the amount of gas mass in the WHIM is 40% of 
total gas averaged over the simulation volumes, in agreement with previous simulations (e.g.. 



Cen & Ostriker 1999 Dave et al. 2001; Cen & Ostriker 2006) and other recent simulations 



(Smith et al. 2011 Dave et al. 2010 Shen & Kelly 2010 Tornatore et al. 2010). For this gas 



we find that the cooling time is tcooi = 6 x 10^ yrs (assuming a temperature of 10^'^K), shorter 
than the Hubble time by a factor > 20. For the strong O VI absorbers considered here, the 
cooling time is still shorter. In other words, either (1) energy is supplied to sustain existing 
O VI gas or (2) new warm gas is accreted or (3) some hotter gas needs to continuously cool 
through the warm phase. Since O VI gas by itself does not define a set of stable systems 
and is spatially well mixed or in close proximity with other phases of gas, this suggests that 
the O VI gas in halos is "transient" in nature. 
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Fig. 4. — shows the metal mass in the warm gas (T = 10^ — lO^K) distributed in the 
density-metaUicity phase space. This is a good proxy for O VI bearing gas. The amount of 
gas mass in the WHIM is 40% of total gas in the simulation volume. 



We consider three sources of warm halo gas: mechanical feedback energy from stellar 
evolution, gravitational binding energy released from halo formation and interactions, and 
direct accretion from the IGM. This simplification sets a framework to make a quantitative 
assessment of these three sources for warm gas that we now describe. We denote Fi, and Fj. 
as the O VI incidence rate (in some convenient units) per blue and red > O.IL^, galaxy due 
to star formation feedback energy heating, Gb and G,- as those due to gravitational heating, 
and Ab and A,, as those due to accreted gas from the IGM. It is useful to stress the distinction 
between G and A. A is gas directly accreted from IGM that is either already warm or heated 
up to be warm by compression upon accretion onto the halo. On the other hand, G is gas 
that is shock heated to the warm phase or to a hotter phase that cools back down to become 
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Fig. 5. — shows the cumulative gas mass as a function of radius for cold (dashed curves), 
warm-hot (solid curves) and hot gas (dotted curves) around blue (blue curves) and red 
(red curves) galaxies at z = 0.2. Also shown the the black triangle is the lower limit from 



observations of Tumlinson et al. (2011a) for star forming galaxies. The horizontal blue 



and red dot-dashed lines are the amount of warm gas the respective star formation rate 
can possibly produce. Additional data point for cold (T < lO^K) gas in early-type galaxies 
within ISOkpc is also plotted as the the red square with the errorbars indicating an estimated 



vertical range from observations of Thom et al. (2012) based on 15 early-type galaxies. 



warm. Restricting our analysis to within a galactocentric radius of ISOkpc and reading off 
numbers from the red curve in Figure |3| we obtain two relations: 



Fk + Gb + Ah 



F, + G,. + A, = 1. 



(1) 



An additional reasonable assumption is now made: feedback heating rate S,. (Sb) is propor- 
tional to average star formation rate SFRj (SFRb), which in turn is proportional to their 
respective gas accretion rate A,. (Ab). This assumption allows us to lump Fr and A,. {Fb and 
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Sb = Fb + Ak = C^ SFRb 
Sr = Fr + Ar = Cx SFR^, (2) 

where C is a constant. We will return to determine Fr and A,. {Fh and A^) separately later. 
Equation ([I]) is now simplified to: 

Sb + Gb = 5, 

5, + G, = 1. (3) 

The ratio of SFRb to SFR^ can be computed directly in the simulations, found to be 8.4. 
Rounding it down to 8 and combining it with Equation ^ give 

Sb/Sr = 8. (4) 

Lastly, a direct assessment of the relative strength of gravitational heating of warm gas in blue 
and red galaxies is obtained by making the following ansatz: the amount of hot T > lO^K 
gas is proportional to the overall heating rate, to which the gravitational heating rate of 
warm gas is proportional. Figure |5] shows the gas mass of the three halo gas components 
interior to the radius shown in the x-axis. Within the galactocentric radius of ISOkpc it 
is found that the amount of hot halo gas per red > O.IL,, galaxy is twice that of per blue 
> O.IL^, galaxy: 

Gr/Gb = 2. (5) 

Solving Equations ( 3p|[5 ) yields 

5b = 24/5,^5 = 1/5; 
5, = 3/5,^, = 2/5. (6) 

The estimate given in Equations (IS]) is admittedly uncertain. Therefore, an estimate on how 
sensitively conclusions depend on it is instructive. We find that, if we had used Gr/Gb = 1 
(instead of 2), we would have obtained Sb = 32/7, Gb = 3/7,3,- = 4/7,(7,. = 3/7; had we 
used Gr/Gb = 1/2, we would have obtained Sb = 4, Gj, = 1,5", = 1/2, G,. = 1/2. Thus, a 
relatively robust conclusion for the sources of warm halo gas emerges: (1) for red > O.IL* 
galaxies {F,. + A,.) and G,. have the same magnitude, (2) for blue > O.IL* galaxies {Fb + Ab) 
overwhelmingly dominates over Gb- 

It is prudent to have a consistency check for the conclusion that star formation feedback 
may dominate heating of warm gas that produces the observed O VI absorbers in blue galax- 
ies. In Figure [5] the horizontal blue dot-dashed line is obtained by assuming a Chabrier-like 
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IMF that is used in the simulations, which translates to 2/3 x lO^^SFR x tcooi x c^ / [klQ^-^ K) , 
where c is speed of light and k Boltzmann constant; also assumed is that 2/3 of the initial 
supernova energy is converted to gas thermal energy, which is the asymptotic value for Se- 
dov explosions, SFR the respective average star formation rate per > O.IL^, blue galaxy, 
i^cooi = 6 X lO^yrs an estimated cooling time for warm halo gas. From this illustration we 
see that with about 20% efficiency of heating warm gas, star formation feedback energy is 
already adequate for accounting for all the observed warm gas around blue galaxies. We 
therefore conclude that the required energy from star formation feedback to heat up the 
warm gas is available and our conclusions are self-consistent, even if the direct accretion 
contribution is zero, which we will show is not. Our results on warm gas mass are also in 



reasonable agreement with observations of Tumlinson et al. (2011a), so is the oxygen mass 



contained in the warm component, as shown in Figure [6] 

Let us now determine E,. and A,- {Fb and A^) individually in the following way. We 
compute warm metal mass that have inward and outward radial velocities within a radial 
shell at r = [50, 150]kpc separately for all red > O.IL* galaxies and all blue > O.IL,, galaxies, 
denoting inflow warm metal mass as Mz{vr < 0) and outflow warm metal mass as Mz{vr > 
0), where Vr is radial velocity of a gas element with positive being outflowing and negative 
being inflowing. We define the infiow warm metal fraction as /in = Mz{vr < 0)/{Mz{vr < 
0) + Mz{vr > 0)), which is listed as the first of the three elements in each entry in Table 
1 under the column \vr\ > Okm/s. We also compute the mean metallicities (in solar units) 
for the infiow and outfiow warm gas, which are the second and third of the three elements 
in each entry in Table 1. Four separated cases are given: (1) red galaxies in C run (C red), 
(2) red galaxies in V run (V red), (3) blue galaxies in C run (C blue), (4) blue galaxies in V 
run (V blue). In order to make sure that infiow and outfiow are not confused with random 
motions of gas in a Maxwellian like distribution, we separately limit the magnitude of infall 
and outfiow radial velocities to greater than 100 km/s and 250 km/s, and listed the computed 
quantities under the third column \vr\ > 100 km/s and the fourth column \vr\ > 250 km/s, 
respectively. 

Table 1. Warm Infiow and Outfiow at r = [50 — 150] kpc Radial Shell 





f r > km/s 


\vr\ > 100 km/s 


Vr > 250 km/s 




(jim ^in/ ■^0; ^out/ ■^o) 


(/in; ^in/ ■^0; ^out/ ■^o) 


(jim ^in/ ■^0; ^out/ ^q) 


Cred 


(58%, 0.27, 0.17) 


(58%, 0.29, 0.17) 


(59%, 0.31, 0.17) 


Vred 


(51%, 0.21, 0.29) 


(61%, 0.18, 0.26) 


(65%, 0.11, 0.33) 


Cblue 


(54%, 0.099, 0.10) 


(55%, 0.099, 0.10) 


(55%, 0.099, 0.10) 


Vblue 


(52%, 0.10, 0.14) 


(52%, 0.09, 0.16) 


(46%, 0.08, 0.24) 
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Fig. 6. — shows the cumulative metal mass as a function of radius for cold (dashed curves), 
warm-hot (solid curves) and hot gas (dotted curves) around blue (blue curves) and red (red 
curves) galaxies at z = 0.2. 



It is interesting to first take a closer look at the difference in metallicities between 
inflow and outflow gas. The warm inflow gas in red galaxies in the C run has consistently 
higher metallicity than warm outflow gas, Zm = (0.27 — 0.31) Z0 versus Zout = 0.17 Z0. 
The opposite holds for red galaxies in the V run: Zin = (0.11 — 0.21) Z© versus Zout = 
(0.26 — 0.33) Z0. The warm inflow gas in blue galaxies in the C run has about the same 
metallicity as warm outflow gas at Z = (0.09 — 0.1) Zq. The warm inflow gas in blue 
galaxies in the V run, on the other hand, has a substantially lower metallicity than the warm 
outflow gas, Zin = (0.08 — 0.10) Z© versus Zout = (0.14 — 0.24) Zq. Except in the case of C 
blue, we note that the inflow and outflow gas has different metallicities, with the difference 
being larger when a higher flow velocitiy is imposed in the selection. This difference in 
metallicity demonstrates that the warm inflows and outflows are distinct dynamical entities, 
not random motions in a well-mixed gas, making our distinction of inflows and outflows 
physically meaningful. A physical explanation for the metallicity trends found can be made 
as follows. In low density environment (i.e., in the V run) circumgalactic medium has not 
been enriched to a high level and hot gas is not prevalent. As a result, warm (and possibly 
cold) inflows of relatively low metallicities still exist at low redshift. The progression from 
blue to red galaxies in the V run reflects a progression from very low density regions (i.e., 
true voids) to dense filaments and group environments, with higher metallicities for both 
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inflows and outflows in the denser environments in the V run; but the difference between 
inflow and outflow metalhcities remains. For red galaxies in high density environments (C 
run) the circumgalactic medium has been enriched to higher metalhcities. Higher cooling 
rates of higher-metallicity gas in relatively hot environments preferentially produces higher- 
metallicity warm gas that originates from hot gas and has now cooled to become warm gas. 
The blue galaxies in the C run are primarily in cosmic fllaments and the metallicity of the 
inflow gas is about 0.1 Z©, which happens to coincide with the metallicity of the outflow 
gas. One needs to realize that at the radial shell r = [50 — 150]kpc over which the tabulated 
quantities are computed, the outflow gas originated in star forming regions has loaded a 
substantial amount of interstellar and circumgalactic medium in the propagation process. 

Let us now turn to the warm inflow and outflow metal mass. It appears that the 
fraction of inflow warm metals (out of all warm metals) lies in a relatively narrow range 
/in = 45 — 65%. For our present purpose we will just say F,. = A,- and Fi, = Ab. Armed with 
these two relations our best estimates for various contributions to the observed warm halo 
metals, as a good proxy for the O VI absorption, can be summarized as follows. 

• For red > O.IL,,, galaxies at z = 0.2 contributions to warm metals in the halo gas from 
star formation feedback (F,.), accretion of intergalactic medium [A,.) and gravitational shock 
heating (G,) are (F,,A,,a) = (30%, 30%, 40%). 

• For blue > O.IL^, galaxies ai z = 0.2 contributions to warm metals in the halo gas from 
the three sources are {Fb,Ab,Gb) = (48%, 48%, 4%). 

• Dependencies of warm halo gas metalhcities on galaxy type and environment are com- 
plex but physically understandable. For red galaxies, the metallicity of inflowing warm 
gas increases with increasing environmental overdensity, whereas that of outflowing warm 
gas decreases with increasing environmental overdensity. For blue galaxies, the metallicity 
of inflowing warm gas depends very weakly on environmental overdensity, whereas that of 
outflowing warm gas decreases with increasing environmental overdensity. As a whole, the 
mean metallicity of warm halo gas in red galaxies is ^ 0.25 Zq, while that of blue galaxies is 
^O.IIZq. 

We suggest that these estimates of source fractions are not seriously in error on average, 
if one is satisfled with an accuracy of a factor of two. The relative metallicity estimates 
should be quite robust with errors much smaller than a factor of two. It is stressed that 
these estimates are averaged over many red and blue galaxies and one is not expected to 
have been led to think that the correlations (such as between warm gas mass and SFR) hold 
strictly for individual galaxies. Rather, we expect large variations from galaxy to galaxy, 
even at a flxed star formation rate. Figure [7] makes this important point clear, which shows 
that, while there is a positive correlation between warm metal mass within 150kpc radius 
and SFR for galaxies with non-negligible SFR (i.e., appearing in the SFR range shown), a 



dispersion of ~ 1 dex in warm metal mass at a fixed SFR in the range of 0.1 — 100 Mq yr^^ 
exists. The goodness of the fit can be used as a way to rephrase this significant dispersion. 
If one assumes that the errorbar size is each log mass determination for each shown galaxy 
is 1, one finds that the chi-square per degree of the fitting line (green) is 0.80, indicating 
that the correlation between logMz(T = 10^~^K) and log SFR is only good to about 1 dex 
in warm metal gas mass. 




log SFR (Msun/yr) 

Fig. 7. — shows the metal mass in the warm gas Mz within a galactocentric radius of 
ISOkpc as a function of the SFR of the galaxy at ^ = 0.2. Each red dot is a galaxy. The 
green curve shows the best linear regression, logMz(T = 1O^~^K)/M0 = 0.32 log SFR + 7.1, 
for the galaxies shown. 



3.4. Composition of Low-z Halo Gas 



In Cen (2012a) we show that the properties of O VI absorption lines at low redshift. 



including their abundance, Doppler-column density distribution, temperature range, metal- 
licity and coincidence between O VII and O VI lines, are all in good agreement with obser- 



vations (Danforth & ShuU 2008 Tripp et al. 2008 Yao et al. 2009). In the above we have 



shown that O Vl-galaxies relations as well as oxygen mass in galaxies in the simulations are 
also in excellent agreement with observations. These tests together are non-trivial and lend 
us significant confidence to now examine the overall composition of halo gas at low-z. 



Figure p^ shows the differential (left panel) and cumulative (right panel) mass fractions 
of each gas component as a function of galactocentric distance for red (red curves) and blue 
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Fig. 8. — shows the differential (left panel) and cumulative (right panel) gas mass fractions 
as a function of radius for cold (dashed curves), warm- hot (solid curves) and hot gas (dotted 
curves) around blue (blue curves) and red (red curves) > 0.1^* galaxies aX z = 0.2. 
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Fig. 9. — shows the differential (left panel) and cumulative (right panel) gas metals mass 
fractions as a function of radius for cold (dashed curves), warm-hot (solid curves) and hot 
gas (dotted curves) around blue (blue curves) and red (red curves) galaxies at z = 0.2. 



(blue curves) galaxies. We note that the fluctuating behaviors (mostly in the differential 
functions on the left panel) are due to occasional dense cold clumps in neighboring galaxies. 
Overall, we see that within about (10,30)kpc for (red, blue) galaxies the cold (T < lO^K) 
gas component completely dominates, making up about (80%, > 95%) of all gas at these 
radii. For both > O.IL^, (red, blue) galaxies cold gas remains the major component up to 
r = (30, 150)kpc, within which its mass comprises 50% of all gas. At r > (30, 200)kpc for 
(red, blue) galaxies the hot (T > lO^K) gas component dominates. The warm gas component. 
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while having been extensively probed observationally, appears to be a minority in both red 
and blue galaxies at all radii. The warm component's contribution to the overall gas content 
reaches its peak value of ~ 30% at r = 100 — 300kpc for blue galaxies, whereas in red galaxies 
it is negligible at r < lOkpc and hovers around 5% level at r = 30 — lOOOkpc. The prevalence 
of cold gas at small radii in red (i.e., low star formation activities) galaxies is intriguing and 
perhaps surprising to some extent. Some recent observations indicate that early- type galaxies 
in the real universe do appear to contain a substantial amount of cold gas, consistent with 



our findings. For example, Thom et al. (2012) infer a mean mass of 10^ — 10^^ M0 of gas 
with T < lO^K at r < ISOkpc based on a sample of 15 early-type galaxies at low redshift 
from COS observations, which is shown as the red square (its horizontal position is slightly 
shifted to the right for display clarity) in Figure [SJ Their inferred range is in fact consistent 
with our computed value of ~ 6 x 10^° Mq of cold T < lO^K gas for red > O.IL,, galaxies 
shown as the red dashed curve in Figure [5j 

Figure [9] that is analogous to Figure [8] shows the corresponding distributions for metals 
mass fractions in the three components for red and blue galaxies. The overall trends are 
similar to those for total warm gas mass. We note one significant difference here. The 
overall dominance of metals mass in cold gas extends further out radially for both red and 
blue galaxies, whereas the contributions to metal mass from the other two components are 
compensatorily reduced. For example, we find that the radius within which the cold mass 
component makes up 50% of total gas in (red, blue) galaxies is (40,150)kpc, where the 
radius within which the cold mass component makes up 50% of total gas metals in (red, 
blue) galaxies becomes (200,500)kpc. This is largely due to a significantly higher metallicity 
of the cold component in both red and blue galaxies compared to the other two components. 



as shown in Figure 10 We also note from Figure 10 that the warm gas in red galaxies 
has a higher metallicity than in blue galaxies, as found earlier, with the mean metallicity 
(~ 0.25 Z0, ~ 0.11 Z0) in (red, blue) galaxies within a radius of 150kpc. 
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log r (kpc) 

Fig. 10. — shows as a function of radius the metalhcity for cold (dashed curves), warm-hot 
(sohd curves) and hot gas (dotted curves) around blue (blue curves) and red (red curves) 
galaxies at z = 0.2. 

4. Conclusions 



The distribution and evolution of the intergalactic medium have largely been addressed 
by cosmological simulations (e.g., Cen & Ostriker 1999; Dave et al. 2001 Cen & Chisari 



2011 Cen||2012a ). The global distribution and composition of halo gas in and around galax- 
ies at low redshift [z < 0.5) are addressed here, utilizing state-of-the-art high-resolution 
(460/i~^pc), large-scale cosmological hydrodynamic simulations with validated star forma- 
tion and feedback prescriptions. 

We find that within about (10,30)kpc for (red, blue) > O.IL^, galaxies the cold (T < 
lO^K) gas component is the primary gas component making up about (80%, > 95%) of 
all gas. For both > O.IL* (red, blue) galaxies cold gas remains the major component up 
to r = (30, 150)kpc, within which its mass comprises 50% of all gas. At r > (30,200)kpc 
for (red, blue) galaxies the hot (T > lO^K) gas component dominates. The warm (T = 
10^ — lO^K) gas component makes a minor contribution to the overall gas mass as well as 
gas metal mass in both red and blue galaxies. The warm component's contribution to the 
overall gas content reaches its peak value of ~ 30% at r = 100 — 300kpc for blue galaxies, 
whereas in red galaxies its contribution is negligible at r < lOkpc and capped at 5% level 
at r = 30 — lOOOkpc. Where comparisons with observations are possible, we find that the 
amount of warm gas and and oxygen mass in star forming galaxies are in agreement with 
observations (e.g., Tumlinson et al.|2011b ), so is the amount of cold gas in early- type galaxies 
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at low redshift (e.g., Thorn et al. 2012). The presence of a significant amount cold gas in 
red galaxies at low redshift is new and somewhat surprising. The nature of this cold gas and 
its role in star formation in red galaxies are not addressed in the present paper. This and 
signatures of the predicted dominance of hot gas in blue galaxies (as well as red galaxies) at 
large radii will be addressed elsewhere. 

In addition to the agreement between our simulations and observations with respect to 
the global O VI incidence rate ( Cen||2012a ), we show that our predicted correlations between 



galaxies and strong O VI absorbers of column density A'^ovi ^ 14 cm~^ are in excellent 
agreement with observations with x square per degree of freedom equal to 1.2. On the other 
hand, when we use only photoionized O VI absorbers (T < 3 x lO^K) in the simulations, 
the comparisons between simulations and observations become substantially less favorable 
with X square per degree of freedom equal to 7.6; the significant disagreement stems from 
the photoionized O VI absorbers being too distant from galaxies. The O VI line incidence 
rate per > O.IL,,, galaxy around blue {g — r < 0.6) galaxies in our simulations is higher than 
that around > O.IL^, red {g — r > 0.6) galaxies by a factor of ~ 4 at r < 100 — SOOkpc, 
increasing to > 10 at r < 20kpc, in reasonable agreement with extant observations (e.g.. 



Chen & Mulchaey 2009 Prochaska et al. 2011a), whereas in the photoionization dominated 



model the ratio is zero at r < 100 kpc. Thus, the cross correlations between galaxies and 
O VI lines provide powerful differentiation between collisional and photo ionization models, 
with collisional ionization dominance for strong (iVovi > lO^^cm"^) O VI absorbers being 
favored currently. 

The O Vl-bearing halo gas (i.e., the warm component) is found to be "transient" in 
nature and hence requires constant sources. We perform analysis to unravel the sources 
and their relations to galaxy formation. We find that, on average, to within a factor of 
two, contributions to warm metals in the halo gas from star formation feedback (F,.), accre- 
tion of intergalactic medium [A,.) and gravitational shock heating [G,-) are (i^,, A,., G,.) = 
(30%, 30%, 40%) for red > O.IL* galaxies at z = 0.2. For blue > O.IL, galaxies at z = 0.2 
contributions are (F^, Af,, Gb) = (48%, 48%, 4%). For both red and blue galaxies, the amounts 
of warm gas in inflows and outflows are comparable. The mean metallicity of warm halo gas 
in (red, blue) galaxies is (~ 0.25 Zq, ~ 0.11 Z©). Environmental dependence of O Vl-bearing 
halo gas is as follows. In low density environments the metallicity of inflowing warm gas is 
substantially lower than that of outflowing warm gas; the opposite is true in high density 
environments. 
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